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Abstract. - A weakly nonlinear model for two-dimensional Faraday waves over infinite depth 
is derived and studied. Sideband instability of monochromatic standing waves as well as 
non-monochromatic solutions are studied analytically. Persistent irregular regimes are found 
numerically. 

Introduction. - Driven gravity-capillary waves is a well-known pattern-forming system owing 
its popularity to a relative simplicity and diversity of observed phenomena. Despite the wide 
interest in the problem, the theoretical basis of pattern formation in large-aspect-ratio systems 
is not yet fully understood. According to experimental data ^] the complex dynamics are 
observed even at small supercriticalities. In the derivation of a theoretical model one must, 
therefore, consistently incorporate the effect of viscous dissipation, which governs the stability 
threshold and supercritical regimes. Furthermore, description of spatially non-homogeneous 
pattern formation involves "slow" space variables, whose scalings in a weakly supercritical 
domain depend strongly on the way how the dissipation was taken into account. 

Some of earlier models Q ^ describing the envelope dynamics of driven waves in a large- 
aspect-ratio system have been derived from the equations for an inviscid fluid with damping 
incorporated phenomenologically by adding linear damping terms to the amplitude equations. 
As a result, an amplitude saturation due to nonlinear damping, which remains important even 
for low viscosity flj, is thus overlooked. In addition, the model obtained in does not possess 
the symmetries ofthe underlying physical system. The nonlinear damping terms are included 
in some other models j^, R] . However, due to neglecting of the boundary layer in the derivation 
of [^, those terms are of higher asymptotic order than in our derivation (see below). In spite 
of using the consistent multiple-scales method with the viscous boundary layer rigorously 
accounted for, due to the chosen scalings the validity of the model seems to be limited to 
exceedingly small supercriticalities. 

In this letter we derive using results of a single amplitude equation describing the evolution 
of one-dimensional patterns for small supercriticalities. This evolution equation is further 
investigated analytically and numerically. 



Model equations. - In a low-viscosity deep-water large-aspect-ratio system a packet of modes 
with adjacent wavenumbers becomes excited due to the parametric subharmonic resonance and 
locked to the (weak) forcing. In this case the principal part of a 2D motion of the free surface 
can be represented as a superposition of two counter-propagating small and narrow wave packets 

C,{x, t) = ae'^"+'"'* + 6e**— + c.c. (1) 
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Here to is half the driving frequency and the corresponding wavenumber k is obtained from the 
dispersion relation ui^ = gk (E/p)fc^, where E is surface tension, g is the gravity acceleration 
and p is the fluid density. 

In the case of weak forcing and low viscosity the slow spatio-temporal evolution of the envelopes 
a{x,t),b{x,t) is described by the set of equations ^ 

1 ^ - 2 "^2 

9ta - sSi^a = -(7 - -72 )a + i/fc + (ici - pi^)|al a + [ic2 - P2^)\b\ a, (2) 

*"2" * 2" 

dtb + sdxb = -{''i--'y"^)h-ifa-{icx+p-i^)\b\ b - {ic2 + P2^)\a\ b, (3) 

where an asterisk denotes complex conjugate. The coefficients in eqs.(|^,(^) are: a finite group 
speed s = du[k)/dk, a small linear damping 7 and the forcing amplitude /, and real- valued 
nonlinear dispersion coefficients, ci,C2, given by 

1 2uk^ , , EP 3 3(t- 2 3 6a -I- 4 

s — a~\ — , 7 = , J = — 1 — a), a = -, Ci — —a , C2 = — a-| . 

2' ' Lj ' 4^ puj^' 4 3a- 1' 2 3a + 1 

Here / is the amplitude of forcing (in units of g), a describes the nature of the wave: it is a 
gravity wave or a capillary wave in two limiting cases of a = and a — 1, respectively. 

Eqs.(^),(^) were derived in from the Navier-Stokes equations using the multiple-scales 
method. Scalings for amplitudes and their slow variation in space and time were chosen as 
a, 6 ~ 72,9j.a ~ dtci ~ 7. Also, following [^] a stretched coordinate was introduced to resolve 
the boundary layer's dynamics. The time and space variables have been non-dimensionalized 
using u)^^ and k~^ , respectively, as their scales. 

3 

Higher-order corrections to the coefficients of the linear damping terms, 7 2 , and of the cubic 
terms, represent, respectively, dissipation in the viscous boundary layer Q and nonlinear 

interaction between potential wavefield and rotational flow in the boundary layer. Due to their 
cumbersome form, complex-valued coefficients Pi are not presented here explicitly. 

Imaginary parts of the coefficients of the cubic terms are responsible for saturation of the 
amplitude growth by tuning the wave out of the resonance (nonl near frequency shift). Asymp- 
totically small coefficients of the cubic terms are retained in eqs.(2|),(|^ since nonlinear damping 
can become important for small supercriticalities [^, ^, for which a simpler model is derived 
below. 

A model similar to eqs.(|2|),(|^ but without the higher-order damping terms was flrst derived 
using symmetry arguments and used in with the coefficients estimated from the experiments. 
High-order nonlinear terms with real coefficients were included in jsj. However, as a result 
of neglecting of dissipation in the viscous boundary layer, the corresponding coefficients are 
asymptotically smaller (~ 7) than in our model. 

In this paper we consider the case of p = 3fJ(pi +P2) > 0, in which, as we show below, subcritical 
solutions do not exist. We also put aside the case of second-harmonic resonance, a « i, in 
which the coefficients ci,pi diverge indicating a breakdown of the mode. 

Note, that the presence of the first spatial derivatives in eqs.(^, (0) follows only from the scaling 
9a; ~ 9t ~ 7 ~ ,/" which is appropriate for the case under consideration. Indeed, according to 
the linear stability analysis pi, H the neutral curve for / ~ 7 ^ 1 is defined by the equation 

s'^ {1 — k)'^ + = Above the threshold fc ~ 7— §72 the bandwidth of unstable modes grows 
proportionally to / for a wide range of supercriticalities = (/ — fc)/ fc- Hence, the evolution 
of the envelopes takes place on the spatial scale f~^, which makes the spatial derivatives in 
eqs.(|^,(^) of the same asymptotic order as the other terms. A similar "non-traditionally" 
slow spatial variation of amplitude takes place in a parametrically driven oscillatory dissipative 
system 0]. 

We relate the damping parameter and supercriticality as 7 = 7i5^, 71 = 0(1)- The expansions 

dx = S^dx^ -f 5^9.3 + . . . , dt = S^d^ + 5^9*3 + . . . , {a, 6} = {a'^) , fo'^)} + 5i{a^^^ , b^^^} + 
. . . , are then introduced into eqs.(0),(p|) and a hierarchy of problems is obtained and solved 
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at each order. At 0(5^''^) we obtain 

a,, a(i) = a., a'l) = a,, 6(1) = a., 6(1) = 0, a^'^=zb^'K (4) 

At this order we find that the motion has a form of standing waves and the dynamics of their 
envelope evolve on the "slower" time- and space-scales t^jX^ . . .. 

At 0{5^) we obtain using eq.(^) two equations 

_ sd.,a'-''> = -a(2) +i6(2) + ic\a'-'^ \\'-'\ (5) 
a^^fed) +s9^^6(i) = -6(2' -m(2) -»cla(if 6(1', (6) 
where c = Ci + C2. Adding eq.(^) and eq.(^ multiplied by i and using eq.(^) we obtain 

at3a(i' = 9*36(1' =0. 

A solvability condition of eqs.(j^,(^) at order S~ yields an amplitude equation, which in terms 
of the unrescaled amplitude a and space and time variables x,t reads 

^2 isc (? 
dta = (/ - fc)a-p^\a\'^a + —dla H -^dx{\a.fa) H —\a\^dxa - 7r7\o-t"- (7) 

This equation describes the slow evolution of the envelope of 2D small-amplitude driven standing 
waves on the surface of a weakly- viscous fluid. In what follows we consider eq.(|7|) with periodic 
boundary conditions in < x < L. It turns out, that eq.(^ is also relevant for alarge spectrum 
of "^/Z W- limiting cases of very small and "large" supercriticalities there are two 

different mechanisms of amplitude saturation in eq.(n). For supercriticalities (/ — fc)/fc fc 
the viscous cubic term is dominant and the dynamics nas the relaxational form described by the 
Ginzburg-Landau equation |1. In the opposite case of "large" supercriticality (or, equivalently, 
small damping) 7 ^ (/ — ffJ/fc ^ 1 saturation is due to nonlinear frequency shift provided by 
the three last terms in eq.([7|). 

In this last case, eq.(0) becomes non- typical, since the "viscous" cubic term becomes small and 
can be omitted, and after appropriate rescaling there is only one independent parameter. The 
reason for this strong mathematical degeneracy is that the coefficients of the cubic terms in 
eqs.(Q),(^) become purely imaginary, which is a result of a weakly dissipative character of the 
system. In this case the only mechanism of amplitude saturation is the nonlinear frequency shift. 
Along with the linear dispersion (first derivative terms in eqs.(P),(^), it leads to different types 
of bifurcation away from the trivial state to the monochromatic solution a = A{t) exp{i(A: — l)a;)} 
around the resonant wavenumber = 1: supercritical for k > 1, and subcritical for k < 1. This 
implies that the coefficient of the cubic term in the corresponding amplitude equation for A, 
changes its sign at A: = 1, and thus, in the vicinity of the minimum the saturation is provided 
by the quintic term. 

According to eq.(^) the driven standing waves have an abnormally large amplitude, which in the 
typical case of 7 ~ (/ — /c)//c «C 1 is proportional to the quartic root of the supercriticality. This 
peculiar steepness of the off-branching solution constitutes the substantial difference between 
patterns formation in a typical strongly-dissipative system and a parametrically forced weakly- 
dissipative oscillatory system. The second important difference is a non-traditionally slow spatial 
variation of the envelope, dxO, ~ (/ — /c)"^''*, whereas for a dissipative system it is ~ (/ — fc)^^^ 



Analytical and numerical study of the model equation. - We first show that in the subcritical 
domain the equilibrium state is absolutely stable. To prove this we multiply eq.(^ by a* and 
integrate over [0, L]. Adding then the complex conjugate of the latter and integrating by parts 
we obtain 

^dt J \a\^dx^ j^^ (^f - fc)\a\^ - p^J^\a\* - ^\isdxa ~ c\a\'^a\^^ dx 



that yields a decay of perturbations when f < fc and p > 0. 
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Assuming now p > 0, / > fc we rescale eq.(|7|) using new variables r = tp^ /2c? , ^ = 
xpf ^ / sc, a = a\c\/ yjpf In a rescaled form eq.(0) reads 

drd — ea-\- d^a — 2\a\'^a + 2il3d^{\a\^)a + 2i\a\^d(a — \a\'^a. (8) 

Eq.(^ contains only two parameters: the efltective supercriticality e — [(/ — fc)/fc]{2(?/'yp^), 
and /3 = ci /c, whose value depends on a and varies in a wide range. 

We now consider the simplest non-trivial solution of the form a(^, r) = A{t) exp(iif^), describ- 
ing a standing monochromatic (MC) wave with the modified wavenumber fc = 1 -)- Kpf ^ / sc. 
In this case the time evolution of the amplitude has the gradient form 

drA = -dAF, where F=^{K + A^f + ^A*^^-A\ (9) 

Since F is bounded from below for fixed K, any MC perturbation results either in the trivial 
solution a = or in one of the stationary solutions 

a±=A±e'^^, A% = -K -1±V£ + 2K + 1. (10) 

It follows from eqs.(|^,(^o|) that the solution a_, which exists for K < —1 and branches off the 
trivial solution subcritically, is unstable. 

The bifurcation diagram of standing MC wave is shown in fig. la for different e. It is seen that 
for large values of e, such that ef <^ 1 (this condition supports the asymptotic method), a 
solution with an arbitrarily large amplitude is possible. If such a solution materializes as a 
result of the time evolution of a finite perturbation \a\ ~ 1, then this would mean the failure of 
the model eq.(^. However, as our numerical study shows (see below), there is a mechanism of 
the amplitude saturation owing to the subcritical character of bifurcation of the large-amplitude 
solution a+ for K < —1. 

We now proceed to sideband (SB) stability analysis of the MC solution a+. The perturbed 
solution is expressed as 

a = e*^« {A+ + 016'"^+^" -I- 026-'"^+^*"). (11) 
growth rate. Upon substitution of eq.(^) into eq.(^) we find 
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A = -q^B{B + K + 1)± ^^B^iB + 1 + + q'^{K + B){K + B + 2/3B), where B = A 

The instability is long-wave and the threshold is defined by the minimal wavenumber allowed 
by the periodic constraint. For an infinite system qmin — *■ 0, and the threshold is defined by 
the equation B = {B + K){K + 2/35). For positive /3's the SB-stability domain is adjacent 
to the saddle- node bifurcation curve esiv = —2K — 1 (curve 1 in fig. lb), and shrinks when /3 
increases. For negative /3's the SB-stability domain widens and moves towards the center of the 
linear instability domain bounded by the neutral curve s = (dotted curve in fig. lb). 

We now show an existence of a stationary non-MC solution of the form a — A{Q exp{i9{^)). 
For this we separate eq.(|^ into real and imaginary parts to obtain 

eA + d^A - 2A^ - {dfOfA - 2A^d(_d - = 0, (12) 
2di,Ad(_e -f- Adle + 2(2/3 -I- \)A'd^A = 0. (13) 
Eq.(p^) can be integrated allowing an elimination of the phase Q from eq.(^2|): 

a|^ = -a,4t/, whereC/=^-t-^(£ + 2C/3-C)-iA*-^(l-2/3)^ (14) 

This equation describes a particle motion in a potential field, where ^ is interpreted as time. 
The potential \J has a local minimum at ^ = . A periodic motion of the particle around the 
minimum corresponds to a non-MC solution periodic in space. Two integration constants can 
be determined from the periodicity conditions for the amplitude A and the local wavenumber 
d^Q. In the presence of the periodic constraint one finds a countable family of periodic solutions, 
which are spatially non-MC. In an unbounded system (in the absence of periodic constraint) the 
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continuum of solutions exists, including a stationary solitary wave corresponding to a homoclinic 
orbit in the phase plane of eq.(Q. 

We have also investigated numerically Eq.(^) along with periodic boundary conditions to study 
more complex non-stationary regimes. We used a pseudo-spectral method with a second-order 
time-marching and up to 256 spatial Fourier modes. In a full agreement with our stability 
analysis, the SB-unstable MC solution breaks down if it is initially perturbed by a sideband 
perturbation. The resulting regime depends on whether for a given e, all SB-stable solutions are 
deeply subcritical {e.g. /? = 2, e = 8 in fig. lb) or there exist SB-stable solutions in a supercritical 
domain {e.g. /? = — 1 in fig. lb). 

In the first case the breakdown of an unstable solution leads (via a phase-slip in the underlying 
standing waves) to onset of the stationary MC waves with a diff'erent wavenumber from the 
SB-stability domain (see fig. 2a). This is a typical behavior for a negative or slightly positive /3, 
for which the domains of SB-stability and instability of the trivial solution overlap for moderate 
e. Non-MC solutions studied above analytically, were not detected in our numerical simulations. 
They seem to belong to unstable manifolds separating the basins of attraction of MC solutions. 

For positive /3 and large values of e, for which SB-stable solutions are in the subcritical domain 
(see the horizontal dashed line in fig. lb), the instability of MC waves typically results in onset 
of an irregular regime. The time-evolution of the Fourier spectrum for 13 — 1,5 = 8 is shown 
in fig. 2b. This irregular behavior was checked using a refined numerical scheme during a rather 
long calculation (r = 10*), in which it remained qualitatively the same. 

In our opinion this irregular behavior is a result of two factors; (i) SB instability of MC solutions 
with moderate K, and (ii) subcritical character of bifurcation of SB-stable MC solutions, which 
does not allow the perturbations with large negative K to attain their stationary values. 

To summarize our results, a new evolution equation describing the weakly supercritical regimes 
of driven surface waves was derived and studied. Its quintic form results in solutions of a 
non-traditionally large amplitude. The sideband instability of a MC wave leads to the onset of 
either a MC wave with a different wavenumber or an irregular motion. 
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Fig. 1. - Bifurcation diagram (a) and stability map (b) of the MC solutions. In (a) stable and 
unstable branches are shown by solid and dashed curves, respectively. The dotted curve in (b) 
denotes the linear stability threshold of the trivial state. The domain of SB-stability of the 
MC solution for /3 = —1 is bounded from below by curve 2, whereas for /3 = 2 it is contained 
between the saddle- node curve 1 and curve 3. 
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Fig. 2. - SB-instability of the resonant {K = 0) standing MC wave leading for /3 = — 1 (a) to 
onset of MC wave with a different wavenumber, and for (3 = 2 (b) to a persistent irregular 
pulsing (shown in Fourier space). The bandwidth of the pulsing is shown in fig. lb by the 
horizontal crossing line. 



